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Ph Abstract: In 2004, two of us proposed a texture, the "Simplest" neutrino mass matrix, 

which predicted sin 6*13 = yj 2Am^ ol / 3 Am^ (m and 5cp = 90°. Using today's measured 

values for neutrino mass-squared differences, this prediction gives sin 2 2#i3 ~ 0.086^o'oo6> 
compared with a measured value of sin 2 2#i3 = 0.093 ± 0.010. Here we present a specific 
model based on S4 symmetry leading to this successful texture in the context of the type-1 
see-saw mechanism, assuming Majorana neutrinos. Similar predictions are obtained relat- 
es) ing #i3 to the light neutrino masses, which are in accord with current experimental limits 
t-h and testable at future experiments. Large CP asymmetries remain a generic prediction of 
the texture. 
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1 Introduction 

Leptonic mixing is characterised by two large mixing angles, Q\i ~ 35° and #23 — 45°, and 
one small angle, #13. For several years, the data on neutrino oscillations were compatible 
with #13 = 0, and the data together were approximated by the tribimaximal (TBM) mixing 
matrix, proposed in 2002 [1]. TBM has been used by many authors as a starting point for 
model building. The Daya Bay Reactor Neutrino Experiment [2] has recently measured 
the value for the reactor mixing angle, sin 2 2^13 = 0.092 ± 0.016 (stat.) ± 0.005 (syst.), 
showing definitively that #13 is non-zero. Similarly, the RENO Experiment [3] also made 
a compatible measurement, sin 2 2#i3 = 0.113 ± 0.013 (stat.) ± 0.019 (syst.). In models of 
leptonic mixing based on discrete symmetries starting with TBM, it is possible to generate 
non-zero #13 by introducing higher order corrections, but in the generic case, the deviations 
produced should be of the same order for all three mixing angles [4]. However, since the 
experimentally allowed deviation of #12 from its TBM value, sin 2 #12 = 1/3, is small, it is 
difficult to generate the relatively large experimental value of #13 in this way. 

Anticipating an eventual non-zero value for #13, two of us proposed several general- 
isations of the TBM texture in 2002 [5], in which the condition #13 = was relaxed in 
various ways. For example, in a basis in which the charged lepton mass matrix is diagonal, 
a hermitian neutrino mass matrix which has /U-r reflection symmetry [6] and democracy 
(each of its rows and columns sums to a common value [7]) leads to "trixmaximal" mixing 
(T%M) [5]. The most general hermitian mass matrix giving T%M takes the form: 
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where k, a, b, and c are real parameters. 
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The mass matrix arising from the Majorana mass term for the neutrinos should be 
complex-symmetric. Therefore we wish to determine the general complex-symmetric mass 
matrix that generates T^M. One way to achieve this is to multiply the matrix Mjj by the 
li-t exchange operator P to get the complex-symmetric matrix M$: 
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where P = ( o o l 
V o l o , 

The only free mixing angle in T^M is the angle #13 which is uniquely determined by 
the real parameter k in the matrix Ms- The four mixing observables in this texture are 
given by: 

|t/d 2 4 **h.-l, sin^ 13 = 1(1.-4=); S0P= i' (1 ' 3) 

We remark that with 5cp = 90°, T%M corresponds to maximal CP violation, for given 
values of the three real mixing angles. The T^M mixing matrix, U, defined by Eq. (1.3), 
which diagonalises the matrix Ms of Eq. (1.2), also diagonalises each of its three terms 
independently. The first term of Ms alone is in fact, sufficient to generate T%M, having 
three distinct eigenvalues. The second and the third terms of Ms give two degenerate and 
three degenerate eigenvalues respectively. As with any complex-symmetric matrix, Ms can 
be diagonalised using a unitary matrix and its transpose to give real positive eigenvalues. 
In other words: 

U ] M S U* = DiagflaVl + 3k 2 - b + c\,\a + 2b + c\,\ - a^/l + 3k 2 - b + c|). (1.4) 

A special case of the T%M texture, known as "Simplest" neutrino mixing was proposed 
in 2004 [8] (after having been introduced and discussed briefly already in 2002 [5]) by setting 
b = in the texture of Eq. (1.1). Its eigenvalues are given by the RHS of Eq. (1.4) with 
6 = 0. Simplest neutrino mixing yields an exact and very straightforward relation between 
the reactor mixing angle (see Eq. (1.3)) and the eigenvalues, ef. 

sin^ *£^4 (1.5) 
3 (e 3 - ei) 

In the original publications [5, 8], this texture was proposed for M 2 := M^M^, in which case 
the eigenvalues are the neutrino masses-squared, resulting in the very successful prediction: 



^ = \IP^- (1-6) 

V 3 Am atm 



ie. sin 2 2(9i 3 = 0.086±g;gg| (Predicted in 2002/2004 [5, 8]) (1.7) 
cf. sin 2 2(9i 3 = 0.093 ± 0.010 (Measured in 2012 [2, 3]). (1.8) 
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In view of the very encouraging phenomenological success of the "Simplest" mixing 
hypothesis, we here propose a model for it based on the symmetric group of degree four 
(S4). In the model, however, we make the following two changes with respect to the origi- 
nal hypothesis: the "Simplest" mass matrix form is adopted for the mass matrix itself (as 
opposed to its hermitian square); in order to exploit the type-I see-saw mechanism a Ma- 
jorana mass term is assumed (coupling between two heavy right-handed neutrinos). Thus, 
we get a Majorana neutrino mass matrix of the following "Simplest" complex-symmetric 
form: 

/ 1 -ik ik \ (l ()\ 
M„(Majorana) = a -ik 1 + ik + c 1 . (1.9) 
yik l-ik) \0 1 Oy 

In the following section we construct such a Majorana mass matrix assuming symmetry 
under the S4 group. The neutrino Dirac mass matrix (coupling between the left-handed 
and the right-handed neutrinos) is assumed to be proportional to the identity. We show 
that this model has a phenomenology compatible with experiment, and we use it to predict 
the masses of the light neutrinos. 



2 The group S4 and the [i-t rotated basis 

54, the group of permutations of four objects, is the symmetry group of the cube and the 
octahedron. In an abstract form, the group has the presentation 

{a,b\a 2 = b 3 = {abf = e), (2.1) 

where a, b and ab represent the orientation-preserving rotations of a cube through angles 
7r, ^ and I respectively. This is shown in Fig. 1, where axis a and axis?, correspond to the 
generators a and b respectively. The character table for the S4 group can be found in [9]. 
We may work in a basis with: 
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as is used frequently in model building. In this basis the coordinate system is oriented such 
that the coordinate axes are normal to the faces of the cube, as shown in Fig. 1(a). So the 
x, y and z coordinate axes are the symmetry axes of ^-rotations. In a model constructed 
with a neutrino triplet (z/ e , z/^, z/ r ) defined parallel to the coordinate axes (x,y,z) in the 
above basis, v e , and v T are simply the invariant eigenstates (eigenstates with eigenvalue 
equal to +1) of these particular ^-rotations . This choice of eigenstates is straightforward, 
and is the one used in most models using this group so far. It is however, by no means 
the only choice, and there is no reason why we should not define the flavour basis states 
in a different way. To construct a model for deriving the "Simplest" texture, it will prove 

1 Abstractly they are also the eigenvectors of the corresponding elements of the three-element conjugacy 
class C3 (-^-rotations about axes passing through face centres). 
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(a) (b) 
Figure 1. Octahedral symmetry in the /i-r rotated basis (x',y',z'). 



useful to define the and v r basis states rotated by an angle \ relative to the x, y and z 
coordinate axes defined above, using the rotation matrix, 



R 




(2.3) 



The matrix R represents a rotation about the x axis by an angle 7r/4, relative to the cube. 
We also rotate the y and z coordinate axes to align with the new and v T flavour basis 
states respectively. In the rotated coordinate system, (x',y',z') in Fig. 1(b), we have the 
group generators: 



a 



R.a.R? 




b -> R.b.R^ 




(2.4) 



The state v e is unchanged and still corresponds to the ^-rotation symmetry of the cube 
about the x-axis. However, and v T are no longer invariant eigenstates of ^-rotations, 
but rather are invariant under n rotations of the cube, as may be seen in Fig. 1(b). We 
call this new basis, the \x-r rotated basis. 

The three-dimensional representation of 5*4 corresponding to the rotational symmetries 
of the cube is denoted by 3' here. Thus the neutrino triplet (u e , u^, u T ) belongs to the 3' 
representation. A Majorana mass term contains two neutrino fields, and thus it is of 
interest to consider the tensor product decomposition of two 3's. This decomposition is as 
follows: 



3'x3' = l + 2 + 3 + 3', 



(2.5) 
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where 1 is the trivial representation, and in the fi-r rotated basis we have: 



Xi 



X3 



— (v e .v e + v^.v^ + v T .v T ), 
-±=(l/ T .l/ e + v e .u T ) I , 



X2 



X 3 



±JyV T .V e - V e .V T ) 



0, 



(2.6) 



where the bi-linears xi, X2, X3, X3 transform as 1, 2, 3, 3' respectively. The product Vi.Vj 
is the Lorentz invariant product of the right-handed neutrino Weyl spinors. Obviously the 
terms in x'3 m Eq. (2.6) vanish. 

We now assume three types of flavons, <pi, cj>2 and <^3 which transform as 1, 2 and 3 
respectively, allowing us to write an invariant mass term: 



Inv 



cixi^i + C2X2V2 + C3X3 h, 



(2.7) 



where c\, C2 and C3 are constants. Once the flavons acquire specific forms of vacuum 
expectation values (vevs), the required mass matrix can be obtained from the invariant 
mass term given in Eq. (2.7). Suppose the flavons get vevs 
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the mass matrix obtained will be 



M = c\I + 
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(2.9) 



■1. 



which is in the same form as Eq. (1.9), assuming c\ and C2 to be real and C3 to be imaginary 
(where here the trace is removed from the second term and is absorbed into the first). The 
details of how to obtain the above mentioned flavon vevs using the minimisation of flavon 
potentials are given in Appendix A. 



3 The model 



The model is constructed in a Standard Model framework with the addition of heavy 
right-handed neutrinos. Through a type-1 seesaw mechanism, light Majorana neutrinos 
are produced. The fermion and flavon content of the model with representations to which 
they belong is given in the Table 1. We construct a diagonal mass matrix for the charged 
leptons using the mass term given by Eq. (3.1). The C2 flavour symmetries Cie-, C^l and C2 T 
ensure that the flavons 4>' 3e , <p' 3fl and 4>' 3t couple to only en, fiR and tr respectively. In order 
to construct the Dirac mass term containing the right-handed neutrinos and the left-handed 
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lepton doublets, we postulate a singlet flavon 07 . We also introduce another C2 symmetry 
C2D to allow only the singlet flavon 07 to enter the Dirac mass term, Eq. (3.2), and thus 
to make the neutrino Dirac mass matrix proportional to the identity. The Majorana mass 
term containing the right-handed neutrinos, Eq. (3.3), leads to a mass matrix of the form 
given by Eq. (2.9) as explained in the previous section. The Standard Model Higgs field is 
assigned to the trivial representation of S4. 
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Table 1. The flavour structure of the model. L are the three left-handed lepton weak isospin 
doublets and vr are the three right-handed heavy neutrinos. The 3' representations are in the 
fi-r rotated basis and the 2 and the 3 representations are in the basis given by the tensor product 
expansion in Eqs. (2.5,2.6). 



For the charged leptons, the mass term is of the form 



y e tf e R (j)' Ze + y^tf fi R cp' 3fM + j/ r L t r R / 3r J ^ + H.C. 



H 



(3.1) 



where H is the standard model Higgs, A is the cut-off scale and the yi are coupling con- 
stants. After the spontaneous breaking of the weak gauge symmetry and the flavour sym- 
metry with the flavons, 4>' 3e , 0' 3m and 3r , getting vevs of (1,0, 0) T , (0, 1, 0) T and (0, 0, 1) T 
respectively, we get the required masses m e , and m T for the charged leptons. The Dirac 
mass term for the neutrinos takes the form 



y w LK R ^j-H + H.C. 



(3.2) 



where H is the conjugate Higgs and y w is a coupling. We also have the Majorana mass 
term for the neutrinos: 

,Tj._ , ,.._..T^ 1 (33) 



(?/lXl01 + 2/2X2 02 + iy3X3 03) V, 

where the Xi are given by the expressions in Eqs. (2.6) and the yi (i = 1, 2, 3) are couplings 
leading to very heavy right-handed Majorana masses. The flavons 0i, 02 and 3 getting 
the vevs given in Eqs. (2.8) and also having (07) = 1 result in the following 6x6 mass 
matrix M for the neutrinos: 

v T a Mv a (3.4) 



M 



M Dir 




(3.5) 
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with i>i = (v e , v^, v T ) the left-handed neutrino flavour eigenstates where Moi r = j£l and 
Mm & j is of the required form given in Eqs. (1.9) and (2.9). Here, MMaj is at a very high 
mass scale and Moir is of order the weak scale, so that the seesaw mechanism comes into 
play. It can be shown [10] that the effective mass matrix, M ss , generated via the seesaw 
mechanism for the left-handed light neutrinos is in the form 

M ss = -M Dir M^jM Dir . (3.6) 

The matrix ill* diagonalises M ss giving light neutrino masses proportional to ^ 
where e\, ei and e% are the eigenvalues of Mmb.]- 

It is to be noted that we may rotate the fermion and the flavon fields from the fi-T 
rotated basis to the conventional basis and obtain the same physical result, even though 
in the conventional basis the charged-lepton mass matrix will not be diagonal and the 
Majorana mass term for the neutrinos will not be in the "Simplest" texture. Therefore 
the use of \i-t rotated basis in the model ensures that the charged-lepton mass matrix is 
diagonal and the Majorana mass term is in the "Simplest" texture. 



4 Fitting the model with experimental data 

The squared differences of the light neutrinos masses are known experimentally, m\ — m\ = 
75.9 ± 2.1 meV 2 , |m| — m 2 \ = 2430 ± 130 meV 2 . The eigenvalues of the Majorana mass 
matrix, MMaj from Eq. (1.9), are 

ei = c + a\/ 1 + 3k 2 , e2 = c + a, = c — a\/ 1 + 3k 2 (4-1) 

(a, c and k are real in our model, given our earlier assumptions). Since the light neutrino 
masses are inversely proportional to these eigenvalues, we get 



ml- m\ _ 75.9 ±2.1 _ 4 ~ \ 
ml -m 2 2 ~ 2430 ± 130 ~ 



(4.2) 



(r - s) 2 (l- s)(l + 2r + s) 
(r + s) 2 (l + s)(l + 2r-s) 



where r = c/a and s = \/l + 3k 2 . Using Eq. (1.3) we can calculate the parameter k given 
the reactor mixing angle 613. Substituting the value of k in Eq. (4.2) and solving for the 
parameter r, we can predict the values of the light neutrino masses. Eq. (4.2) is cubic in 
r giving three separate real solutions for normal hierarchy and one for inverted hierarchy. 
One of the normal hierarchy solutions gives the wrong sign for the solar mass-squared 
difference, leaving three remaining solutions. Thus the light neutrino masses predicted 
by the model fall into three sets. These results are shown in Fig. 2 where the best fit 
values are used. The error ranges of the mass of the neutrino eigenstate ui for the three 
solutions are shown in Fig. 3. For solution 2 we do not consider a mass above 100 meV in 
order to keep our prediction compatible with the cosmological upper limit of the masses 
of the neutrinos. The WiggleZ Dark Energy Survey [11] gives the strongest cosmological 
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limit so far, m v < 290 meV. It should be emphasised that in trixmaximal mixing 
(T^M), given the three mixing angles #12, #23> 6*13 , maximal CP violation (5cp = 90°) is 
always guaranteed. The positive and negative signs of the parameter k correspond to the 
CP-violating phase 5qp = ±90°. 
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120 - Solution 1 Solution 2 Solution 3 



100 - m 3 

m 2 

80- ~»l 
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40 
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m 2 
mi 

m 2 
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m 3 



Figure 2. The predicted values of the neutrino masses corresponding to the best fit 
(ml - ml = 75.9 meV 2 , |m§ - m\\ = 2430 meV 2 , sin 2 26 13 = 0.098). Case 1 (r = 0.4101) 
and case 2 (r = 14.452) are in normal hierarchy. Case 3 (r = —1.0405) is in inverted hierachy. 




Figure 3. The predicted value of mi (the mass of the neutrino eigenstate v\) vs the measured 
value of sin 2 2#i3. The finite thickness of the bands is due to the errors in the measurement of the 
neutrino mass-squared differences. The red and the black lines indicate the best fit value and the 
errors on sin 2 2# 13 respectively. 
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5 Summary 

We use the "Simplest" texture of neutrino mass matrix to explain the recently measured 
non-zero reactor mixing angle. When used as a hermitian square of the mass matrix, 
this texture successfully predicts the observed reactor angle (with the help of the neutrino 
mass-squared differences). When used as a complex-symmetric Majorana mass term, the 
constraint among the masses and the mixing resulting from the texture can be used to 
predict the unknown lightest neutrino mass (with the help of the reactor angle). We get 
three solutions which are compatible with the measured mass-squared differences. In the 
model we exploit the "p-r rotated basis" , introduced here for the first time, to obtain the 
required texture for the mass matrix. Since the mixing is trixmaximal, the CP-violating 
phase is predicted to be ±f . Large CP violation such as this is potentially testable in 
future experiments. 

This work was supported by the UK Science and Technology Facilities Council (STFC). 
Two of us (PFH and RK) acknowledge the hospitality of the Centre for Fundamental 
Physics (CfFP) at the Rutherford Appleton Laboratory. RK acknowledges support from 
CfFP and the University of Warwick. 



6 Appendix A: Flavon Potentials 

Defining the flavon 02 = (02, 0|), we can construct a third degree invariant, — {4>\) z + 
3(</>2) 2 02- Along with the term (4>\) 2 + (02) 2 (which is basically a U± invariant), we can 
construct the potential 

U(0 2 ) cx -(0*) 3 + 3(^) 2 0i + ((0i) 2 + (^f) 3 (6.1) 

leading to a vev of (fa) = (— \, ^), as required. 

For the flavon 03, it is easier to work in the non-//-r-rotated basis. Defining 3 = 
(03,03,$$) in this basis, £4 invariants can be easily constructed by summing up even 
powers of 03, 2 and 0| symmetrically. A potential constructed in this way is given below: 

U(0 3 ) oc ((0i) 8 + (0|) 8 + (0i) 8 ) + P ((0 3 ) 6 + (0D 6 + + Q ((<Pl) 4 + itif + (4)') 

+ r ((0 3 ) 2 + (01? + (0|) 2 ) + - ((0 3 ) 2 (0|) 2 + {4>l)\4>l? + (6-2) 

for real p, q, r and s under some constraints. This results in (03) = (l,\/2, 0) in the 
non-rotated basis, corresponding to (03) = (1, 1, —1), in the [i-T rotated basis, as required. 
Using the extremisation condition of zero first order derivatives applied at the point 03 = 
(1, y/2, 0), we get the following constraints: 

r = -60 - 21p - 6q, (6.3) 
s = 28 + 9p + 2q. (6.4) 
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To ensure that the extrema points are minima, we need to impose the condition of positive 
definite Hessian matrix. This gives the following inequalities: 



(4 + p) >0, 
(l8 + 5p + q) > 0, 
(208 + 72p + 9p 2 - 8q) > 0. 



(6.5) 
(6.6) 
(6.7) 
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